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A mean-field theory of the spin-Peierls systems based on 
the two dimensional dimerized Heisenberg model is proposed 
by introducing an alternating bond order parameter. Im- 
provements with respect to previous mean-field results are 
found in the one-domensional limit for the ground state and 
the gap energies. In two-dimensions, the analysis of the com- 
petition between the antiferromagnetic long range order and 
the spin-Peierls ordering is given as a function of the coupling 
constants. We show that the lowest energy gap to be observed 
does not have a singlet-triplet character in agreement with the 
low temperature thermodynamic properties of CuGeOa. 

PACS number(s): 75.10.Jm, 75.30.-m, 75.30.Kz, 75.50.Ee 



The recent discovery of a spin-Peierls (SP) transition in 
the inorganic compound CuGe03 prompts renewed inter- 
est for this kind of structural instability.ErH Evidence for 
such a non-magnetic transition has been exemplified in 
several ways. The lattice distorsion has been well estab- 
lished by X-ray and elastic neutron experiments!] - tl The 
magnetic susceptibility decreases exponentially showing 
a gap in the spin excitations .Era This is also confirmed by 
heat capacity measurements which present a thermally, 
activated component below .the SP critical temperature.EI 
As shown by Oseroff et a/.,0 another interesting feature 
is the close proximity between the SP and the antiferro- 
magnetic (AF) ground state for doped CuGeOs samples. 

For the theoretical description of the SP ordered state, 
the alternating Heisenberg-. model has been .extensively 
studied both numericallycm and analyticallyJliJ However 
in the latters, mean-held like decoupling of the Heisen- 
berg interaction in the quasi-fermion representation were 
essentially restricted to the single chain problem. In the 
present work, we introduce a new decoupling for the al- 
ternating bond order parameter which not only improves 
ID mean-field results for the ground state energy and the 
excitation gap but also allows to treat the 2D situation,EII 
namely the effect of interchain Heisenberg exchange in- 
teraction. This turns out to be an important ingredient 
for the interplay between the AF and the SP states. 

We start the analysis with the 2D alternating Heisen- 
berg model 

H = J &2i,j • S2i+l,j + J' /J S2i+1 j • S2i+2j 

J± J]] &i,2j ■ Si^j + l + J± J]] Si^j + l • Si^j + l (1) 
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where J > (J± > 0) and J' > (J' ± > 0) are the 
intrachain (interchain) AF exchange couplings. We then 
follow the mean-field approach given in refs. [12] and 
[13]. In such a treatment, the spin Hamiltonian is trans- 
formed by meaps of the generalized Jordan- Wigner (JW) 
transformation^ 2 ] 
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cjjCij. The Hamiltonian is then written in 



where 

Fourier space by taking into account the bipartite char- 
acter of the lattice, namely 



4 - Mc^c* 



- e (k)cjfc£ + e *(k)cf C £}. 



(3) 



Here A and B label the two sublattices and e(k) = 
j ie ika, _ j 2e - lfe * + J ±lC lk y + J ±2 e- lk y , M = m( J + J' + 
J± + J' ± ) with to, Q, Q 1 , P, P', Ji and J 2 are to be de- 
fined shortly. The Hamiltonian (|J) is diagonalized using 
the following canonical transformation 



Ck = u kfk + Vkdk 



(4) 



where u k = e lak / 2 cos f3k and Vf. = e lQfc / 2 sin/3/j. atf- and 
/3k are given by 



tan at. = 



( J\ - J2) sin k x + (J11 - Jj_2) sin k y 
( Ji + J2) cos k x + (J±i + Jl 2 ) cos k y 



and 



tan(2/? fc )= ±(2Af) -1 x 

( [( Ji + J2) cos k x + ( Jj_i + Jj_ 2 ) cos k y ] 

[( Ji - J 2 ) sin k x + ( J ±1 - J|_ 2 ) sin k y ] 
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where Ji = J(1+2Q), J 2 = J'(1+2Q / ), J±i = J±(l+2P) 
and J±2 = J'±(l + 2P'). As for the dispersion relation, 
it is given by 



£ ± (k) = ±i{Af 2 + | e (k)| 2 } 



1/2 
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where (±) refers to upper and lower band. We have 
introduced the order parameters to = 2(S Z ), Q — 



\( C 2i,j4i+l,j)l Q' = l( C 2mj c L+2j)l' p = \( C u2 3 C 



1 



and P' = |(ci,2j+ict 2j . +2 )| for the staggered magnetiza- 
tion, intra- and interchain bond amplitudes respectively. 

Their equilibrium values are obtained from the mini- 
mization of the total free energy leading to a set of mean 
field equations which can be solved exactly. 

The ground state wave function corresponds to the case 
where the lower band is filled: 



I'iW = IKlO>- 



(6) 



It is formed by the pairs of fermions (c£\ c^) with the 
weights Vk and — u| respectively which correspond from 
(0) to pairs of spins denoted (f, |)k in reciprocal space. 
The ground state is a singlet since ($gs|5' z !^gs} — 
($G5|(S tot ) 2 | < i>Gs) = 0, while excited states with a 
wavevector k', relevant to magnetic susceptibility and 
specific heat etc, consist in creating particle-hole excita- 
tions, namely 

\*ex) = fUv\*Gs) = Y[fld k ,dl\0) (7) 

k 

where the excitation operator can be written as 

/ k dk « u*t£(SS(k) - S z B (V)~(uk) 2 S+(k)S B (k) 

+ (v* k ) 2 S A (k)S+(V, (8) 

where the phase factors in (Q) have been ignored in or- 
der to illustrate that such excitation operator conserves 
Sl ot so that \$ex) is not a triplet. The dispersion re- 
lation for each member of the particle-hole pair is given 
by \{<S>ex\H\<5> E x) - \^gs\U\^gs) = E+(k) with an 
energy gap 



E a = ^{m 2 (J + J' + J x + J' ± ) 2 



\Jl - J2} 2 + [Jxi - J±2} 2 } 1/2 , 



(9) 



at ko = (0, 7r/2). In the case of uniform exchange interac- 
tion, the gap is due to the finite sublattice magnetization 
to, as discussed previously £3 At zero to and J± = J' ± , 
the gap reduces to E^ p — || J\ — J%\ which is due to the 
dimerization along the x-axis. The ground state formed 
by dimers is the so called SP state. An elementary ex- 
citation will correspond to the destruction of a singlet 
(t> l)k from the application of dk - Actually, in one 
dimension one can visualize Eg P for the creation of a 
particle-hole excitation as the energy required to form 
a kink in the dimerized chain. Finally, for longitudinal 
and transverse alternating exchange interactions, the gap 
adds in quadrature and we find for the particle-hole ex- 
citation 



E; 



SP 



;{[Jl- J 2 ] 2 + [J±l- J±2} 2 } 



1/2 
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Static and uniform spin susceptibility x is equivalent to 
the compressibility in the JW representation. Its evalu- 
ation is straightforward and the result is given by 



x(T) = -i,Vi£^ +(k)1 



dE+(k) 



(11) 



where n[x] = (eP x + 1) _1 , g is the Lande factor and 
T>(Eg P ) is the density of states at the energy gap. There- 
fore the lowest energy gap to be observed in experiments 
(e.g. in GuGeOa) like magnetic susceptibilityfiJ and spe- 
cific heatB of the condensed SP state is characterized by 
the above particle-hole (singlet) character. Furthermore, 
in contrast to the critical temperature, the amplitude 
of the zero temperature gap is not predicted to change 
when a magnetic field is applied in agreement with spe- 
cific heatu and acoustic measurementstfl performed under 
low field. 

When both possibilities of SP and AF long range or- 
der are considered (Eq. (^)), the numerical solution of the 
T = mean field equations leads to the phase diagram 
of Fig. [l] for (J'/ J, J±/J) at a fixed J' ± value. Whenever 
to =/= 0, the rotational invariance is broken and the system 
is in the AF state while the line boundary is determined 
when to vanishes. As for the SP phase, it is defined by 
to = and an energy gap solely due to dimerization. In 
the insert of Fig. [j], the magnetization is displayed as a 
function of J'/ J for J± = J' ± = .1J. Therefore when- 
ever the dimerization becomes sufficiently small, the SP 
ordered state is no longer stable and a magnetic ordering 
is favored; an increase of the interchain exchanges also 
favors the magnetic ordering. This is consistent with 
the situation found in real systems like quasi-lD organic 
materialsEj where the application of hydrostatic pressure 
is well known to promote such a crossover. 

Focussing now on the ground state energy Eqs — 
($as\H I^gs) m the to = SP phase, one can compare 
in Figure || the present mean-field approach with the pre- 
vious results obtained by BulaesvskiililJ in the ID limit. 
Thus the choice of an alternated order parameter (Q ^ 
Q') gives rise to a better estimation of the ground state 
energy. In the uniform Heisenberg limit of the model, 
J = J' , the dispersion relation E + (k) = (1 + 2Q)\ sinfc| 
becomes gapless with-J + 2Q w 1.63 in fair agreement 
with the exact result Ji3 

As far as the gap is concerned in this ID limit (insert 
of the Figure 0), the present approach leads to 



E 5 / = E 



C I 1 - J'/ J 



(12) 



with C ~ .8J and the exponent a ~ .71 which is close 
to the Cross and FisheJl3 value (a = 2/3), RG results, 
(a = .76)E3 and exact diagonalization (a = .79 ± .06).E£l 
The present calculation however predicts a finite jump 
Eq ~ .19 J for the energy gap once | J — J' | is non 
zero; a result not yet confirmed by an exact numerical 
calculation probably due to finite size effects as | J — 

,r h on 

When the effect of interchain exchange coupling is in- 
cluded (Figure |j|) , one can extract in the region where the 
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SP phase is stable (for J±/J < .45, the AF phase is dom- 
inant for all | J — J' |) the universal value a ~ .66 for the 
exponent, while the constant E decreases monotonously 
with J± and C ~ .81J. As a 2D mean- field result, this 
value oLa is closer to the one obtained by Cross and 
FisherB 

In summary, we have proposed a mean-field theory of 
the 2D dimerized Heisenberg model with magnetic and 
alternating bond order parameters for the description of 
the ordered spin-Peierls state and its competition with 
the antiferromagnetic order. In the ID limit, the ground 
state energy and the particle-hole excitation gap profiles 
with dimerization show marked improvements with re- 
spect to previous mean-field results. As far as the exci- 
tation gap is concerned, particle-hole excitations should 
dominate the thermodynamics at low temperature and 
this, consistently with recent measurements on specific 
heat and magnetic susceptibillity for the CuGe03 mate- 
rial. 
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FIG. 1. The phase digram (J'/J,J±/J) is drawn for 
J' ± — .08 J. In the insert, the magnetization m is reported 
as a function of J'/ J for J' ± — J± = .1J. 

FIG. 2. The ground state energy as a function of | J' — J\ 
in the present mean-field approximation compared to the Bu- 
laevskii's results (B). In the insert, the same comparison for 
ID SP gap. 



FIG. 3. The energy gap Eg P as a function of | J — J' | for 
J± = .15 J (curve 1) and ,J± = 0.1J (curve 2). The curve 3 
gives the energy gap which would be obtained in the absence 
of antiferromagnetism (m = 0) for J± — .15 J. In the insert, 
the magnetization is reported for J± — .15 J (curve 1) and 
J± = .1J (curve 2). 



1 M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. Lett. 
TO, 3651 (1993). 

2 M. Hase et al, Phys. Rev .B 48, 9616 (1993) 

3 L. P. Regnault, M. Ain, B Hennion, G. Dhalenne and A. 
Revcolevschi, Int. Conf. Neut. Scat., to appear in Physica 
B. 

4 K. Hirota et al, Phys. Rev. Lett. 73, 736 (1994) 

5 J. P. Pouget et al, Phys. Rev. Lett. 72, 4037 (1994). 

6 T. M. Brill et al, Phys. Rev. Lett. 73, 1545 (1994). 
7 Oseroff et al, Phys. Rev. Lett. 74, 1450 (1995). 

8 W. Duffy and K. P. Barr, Phys. Rev. 165, 647 (1968) 

9 J. C. Bonner and H. W. J. Blote, Phys. Rev. B 25, 6959 
(1982); K. Okamoto, H. Nishimori, and Y. Tagushi, J. 
Phys. Soc. Jp. 55, 1458 (1986); Z. G. Soos, S. Kuwa- 
jima, and J. E. Mihalick, Phys. Rev. B 32, 3124 (1985) 
and references therein; G. Spronken, B. Fourcade, and Y. 
Lepine, Phys. Rev. B 33, 1886 (1986); J. Riera and A. 
Dobry (preprint). 

10 L. N. Bulaevskii, Zh. Eksp. Theor. Fiz. 44, 1008 (1963) 
[Sov. Phys. JETP 17, 684 (1963)]. 

11 N. Katoh and M. Imada, J. Phys. Soc. Jp. 63, 4529 (1994). 

12 M. Azzouz, Phys. Rev. B 48, 6136 (1993). 

13 M. Azzouz, L. Chen, and S. Moukouri, Phys. Rev. B 50, 
6233 (1994). 

14 M. Poirier, A. Revcolevschi, and G. Dhalenne (private com- 
munication) . 

15 See for example, P. Wzietek et al, J. Phys. (Paris) I 3, 171 
(1993) 



3 



